Spinodal decomposition to a lamellar phase: effects of hydrodynamic flow 
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Results are presented for the kinetics of domain growth of a 
two-dimensional fluid quenched from a disordered to a lamel- 
lar phase. At early times when a Lifshitz-Slyozov mechanism 
is operative the growth process proceeds logarithmically in 
time to a frozen state with locked-in defects. However when 
hydrodynamic modes become important, or the fluid is sub- 
jected to shear, the frustration of the system is alleviated and 
the size and orientation of the lamellae attain their equilib- 
rium values. 

PACS numbers: 47.20Hw; 05.70Ln; 83.70Hq 



The kinetics of the growth of ordered phases as a sys- 
tem is quenched from a one- to a two-phase region pro- 
vides a wealth of interesting physical questions [ 0]. Our 
aim here is to present results for phase separation in 
a two-dimensional binary fluid in the case where equi- 
librium corresponds to a lamellar structure where the 
two coexisting phases arrange themselves in stripes. We 
consider both the early-time regime, where the growth 
is driven by a Lifshitz-Slyozov evaporation-condensation 
mechanism [U], and later times (or equivalently lower vis- 
cosities) where hydrodynamic modes alter the kinetic be- 
haviour. Our main conclusions are that when the former 
mechanism is operative the lamellae initially grow loga- 
rithmically but then become frozen in a tangled configu- 
ration. It is only when hydrodynamic movement becomes 
effective that topological defects can be removed allowing 
the lamellae to line up and attain their correct equilib- 
rium width. We show that the frustration of the system 
can also be alleviated by subjecting it to shear. 

Our results are relevant to a wide variety of physi- 
cal systems where competing interactions result in stable 
lamellar phases. One example is a di-block copolymer 
melt comprising A- and B- type chains covalently bonded 
end-to-end in pairs. If the A-B interaction is repulsive 
the tendency towards bulk phase separation at low tem- 
peratures is halted by the covalent bond and a lamellar 
phase can be stabilised [ pi. In microemulsions similar 
structures result from the presence of surfactant at the 
(say) oil-water interface which leads to competition be- 
tween a negative surface tension and positive interface 
curvature free energy [pi. Other examples include dipo- 
lar fluids with long-range Coulombic interactions [|7[] and 
chemically reactive binary mixtures where the ordering 
tendency of the two components is counterbalanced by 
the mixing effect of the reaction [ @]. 

The problem is addressed using lattice Boltzmann sim- 
ulations[S| of the Navier-Stokes and convection-diffusion 
equations for a binary fluid. These simulations follow the 



behaviour of a set of distribution functions, related to the 
fluid densities and velocity, which evolve according to a 
discretised Boltzmann equation. The densities and mo- 
menta are conserved at each step of the simulation thus 
approximating the continuum equations of fluid flow. 

The technique has two advantages which are particu- 
larly relevant here. Firstly it allows simulations on long 
enough time scales to access hydrodynamic behaviour. 
Secondly by using an implementation of the method de- 
scribed by Orlandini et. al. [ |l^ it is possible to introduce 
a free energy into the model such that the fluid relaxes 
to an equilibrium state determined by the minimisation 
of the input free energy. Hence we are able to accurately 
control the equilibrium state that the fluid is trying to 
reach. 

The numerical method is described in detail in refer- 
ence [ noi. Here we consider a model which corresponds 
to an equilibrium free energy density [ O] 
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where ip is & scalar order parameter describing the con- 
centration difference between the two binary compo- 
nents. We take 6, c > to ensure stabihty. For a > the 
fluid is disordered; for a < it prefers an ordered state. 
The nature of this depends on the value of the surface 
tension. For k sufficiently large for fixed c the surface 
tension is positive and the fluid prefers to phase sepa- 
rate into two homogeneous phases. As k, is decreased, 
however, at some n = n^ ^ —0.8 for c ~ \ the surface 
tension becomes negative and it becomes favourable to 
introduce interfaces into the system. The ordered state 
is now a lamellar or striped phase with spacing deter- 
mined by the competition between the free energy gain 
due to interface creation and the repulsive force between 
interfaces. We shall compare the behaviour of the fluid 
under sudden quenches from the disordered to the or- 
dered phases for k greater than and less than Kc- 

Before presenting our results it is useful to summarise 
the known behaviour of a binary fluid quenched into the 
ordered homogeneous phase [ ^. Once domains of the 
two phases are well established experimental and numer- 
ical data indicates that the growth is self-similar with 
a typical length scale R scaling with time as i? ~ t". 
The exponent a is universal depending only on a small 
number of parameters; the dimensionality, conservation 
laws applicable to the order parameter and the presence 
or absence of hydrodynamic degrees of freedom. Three 
growth mechanisms are believed to be operative in two 



dimensions: (i) Lifshitz-Slyozov evaporation of particles 
from interfaces of high curvature to recondense on inter- 
faces of lower curvature with a — 1/3, (ii) a regime with 
a = 1/2 which has been explained as being due to Brow- 
nian diffusion of interfaces and droplets [ E2|, (iii) inertial 
hydrodynamic flow with a = 2/3[ |l3|. Previous lattice 
Boltzmann simulations [ l3 have shown a crossover from 
(i) to (iii) with decreasing viscosity with (ii) not seen as 
might be expected because of the lack of noise in the sim- 
ulations. Lattice gas simulations, which do include noise, 
show a crossover from (ii) to (iii) [ nsl . However puzzles do 
remain. For example, a simulation by Lookman et. al.[ 
M using a Ginzburg-Landau approach, shows a crossover 
from (ii) to (iii) independent of whether or not noise is 
present. 
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FIG. 1. Logarithm of a typical domain size as a function 
of the logarithm of the time for different values of the surface 
tension and a sufficiently high viscosity that the growth pro- 
ceeds by the Lifshitz-Slyozov mechanism. The straight line 
corresponds to a growth exponent a = 1/3. For k sufficiently 
negative the domain size grows logarithmically and then satu- 
rates as the system reaches a frozen state of tangled lamellae. 

Results for the evolution of the domain size R{t) with 
time following a quench to the ordered phase a < for 
different values of n are shown in Fig. 1. R{t) was calcu- 
lated as the first moment of the circularly averaged struc- 
ture factor. These simulations were run with high val- 
ues of the viscosity where the diffusive Lifschitz-Slyozov 
growth mechanism is expected to be operative. 

For K > —0.8, where the surface tension is positive, all 



systems reach the scaling regime (i) after a sufficiently 
long time. As k is decreased the domains take longer 
to reach a given size. This is as expected because the 
driving force for the Lifshitz-Slyozov growth, provided 
by the surface tension, is lower. It is interesting to note 
that the scaling regime sets in at approximately the same 
domain size in each of the runs. The data presented in 
Fig. 1 is from a single run on a lattice of size 128 x 128 as 
our aim here is to provide a qualitative comparison of the 
results as k is decreased rather than the best quantitative 
value for a which has been done elsewhere[ll4J. 

For K < Kc where the system prefers to order in a 
lamellar phase a significantly different behaviour is ap- 
parent. After initial transients there is a region of loga- 
rithmic growth which corresponds to formation of the 
lamellae. The growth then slows down and the fluid 
becomes frozen. A similar pattern of growth has been 
observed for quenches into microemulsion phases [ pT 15 
and for systems with quenched defects [ [l8| . 







FIG. 2. Snapshots of the growth of domains with time. 
Each line, reading from left to right, represents a different 
physical situation: (a) a quench to the homogeneous two 
phase region k > Kc\ (b) a quench to the lamellar phase k < Hc 
in a high viscosity fluid. The lamellae form in a tangled pat- 
tern which becomes frozen in time; (c) a quench to the same 
final parameters as (b) but for a low viscosity fiuid. Hydro- 
dynamic modes allow the lamellae to reorder giving, locally, 
well-defined striped regions. 

A pictorial comparison of the domain growth for k > 
Kc and K < Kc is made in Figs. 2(a) and (b). In the for- 
mer case the system prefers to form elongated domains 
because of the curvature term in the free energy. How- 



ever this does not impede growth into two homogeneous 
phases and the lamellae grow and become more isotropic 
throughout the simulation. For k < Kc, however, (Fig. 
2b), the lamellae form in a tangled pattern which be- 
comes frozen in time. 

We now turn to a primary aim of this Letter which is to 
report the effect of hydrodynamic modes on the quench 
to the lamellar state. To investigate this we ran the sim- 
ulations for K = —0.9 for a low value of the viscosity. 
A pictorial comparison of the domain growth in the long 
time regime for the high and low viscosity cases are given 
in Figs. 2b and 2c respectively. It is apparent that in the 
latter case it becomes easier for the lamellae to grow and 
shrink along their length thus allowing a decrease in the 
frustration in the system. This is because the long-range 
nature of the hydrodynamic modes are allowing the fluid 
to reorder in such a way as to remove topological de- 
fects from the fluid. For the latest accessible times some 
defects remain although the system is locally lamellar. 
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FIG. 3. Typical domain size as a function of the logarithm 
of the time for a quench to the lamellar phase for □ random 
initial conditions and a high viscosity; Q random initial con- 
ditions and a low viscosity; A a sine wave as initial condition 
and a high viscosity. In the second and third cases, for long 
times, the fluid is not frustrated and the lamellae can attain 
their equilibrium width. 

The rearrangements of the lamellae also have the ef- 
fect of allowing them to attain their equilibrium width. 
The growth of the domain size with time following a 
quench to the lamellar phase for high and low viscosi- 



ties is compared in Fig. 3. In the latter case there is 
a relatively rapid decrease in the domain width as hy- 
drodynamic modes become effective. As evidence that 
equilibrium is indeed being achieved we also show in Fig. 
3 a quench for the same (high viscosity) parameters but 
with a one-dimensional sine wave as the initial condition. 
This allows the lamellae to grow easily in an ordered pat- 
tern. The final domain size achieved is very similar to 
that obtained after hydrodynamic rearrangements of the 
frozen state. 

As a further check on the mechanism causing the un- 
locking of the frustrated lamellar state we simulated the 
behaviour of simple topological defects. Examples are 
a broken lamella or an extra portion of lamella inserted 
in an otherwise regular pattern. In each case diffusive 
growth was only able to cause small changes of interface 
shape whereas hydrodynamic rearrangements succeeded 
in removing the defect from the fluid. 




FIG. 4. Configurations of a lamellar fluid when a shear flow 
is applied. The shear velocity is: (a) v=0, (b) v=0.1 and (c) 
v=0.3. 

Finally we consider another mechanism which allows 
topological disorder in the lamellar state to be removed, 
the introduction of a shear. To simulate the shear we 
imposed rigid walls moving with velocity ±v on the left 
and right of the simulation area. A natural way to im- 
plement rigid walls is to use bounce-back conditions at 
the boundary sites [ M. The other edges of the system 
were linked with periodic boundary conditions. 

The fluid was quenched into the lamellar state giv- 
ing the disordered conflguration shown in Fig. 4(a). As 
the shear was applied ordered lamellae quickly formed, 
aligned at an angle to the direction of shear, as shown in 
Figure 4(b). A further increase in the shear rate resulted 
in a decrease in the angle between the direction of the 
lamellae and that of the shear (Figure 4(c)). We believe 
that this angle is determined by a balance between the 
force on the lamellae due to the shear and the elastic force 
which results because they are compressed. Presumably 
on much longer time scales lamellae could be lost from 
the system, releasing the elastic energy and allowing the 
lamellae to turn towards the shear direction. 



Thermal fluctuations are not included in the lat- 
tice Boltzmann simulations and it is important to ask 
whether noise could be responsible for unpinning the 
metastable state before the onset of hydrodynamics. To 
test this we ran a simulation of Langevin dynamics for 
the free energy (P. Tangled lamellar were formed in the 
same way as in the lattice Boltzmann simulations sug- 
gesting that noise is ineffective in bringing the system 
to a global thermodynamic equilibrium. This conclusion 
is supported by cell dynamical calculations [ Q but dis- 
agrees with results based on a time-dependent Ginzburg- 
Landau simulation [p^]. 

To conclude, we have studied domain growth in a fluid 
with a well-characterised, lamellar, equilibrium state. 
Diffusive growth leads to a frozen pattern of tangled 
lamellae. The topological defects in the ffuid can be 
removed by hydrodynamic modes or applied shear. It 
would be interesting to extend this work to three dimen- 
sions where a new physical mechanism, the Rayleigh in- 
stability of tubes, is expected to be operative[|22J. How- 
ever such calculations would be extremely demanding 
computationally. 
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